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SUMMARY 

Thermal neutron flux albedo for concrete due to an incident fast- neutron beam is of 
interest in predicting activation of a nuclear- rated space chamber. Some measurements 
of this quantity were made for limestone concrete by using plutonium-beryllium and 
reactor-beam neutron sources. In both cases, measurements were made by using a 
moving-detector technique to perform the integration of emergent thermal neutron flux. 
The accuracy of results for the reactor source was limited by lack of beam- spectrum 
data, but a value of flux albedo for the PuBe source was determined to be 0. 13±0. 04. 


INTRODUCTION 

Neutron albedo from thick concrete slabs is of concern in testing of nuclear-powered 
vehicles in a large space environment chamber now under construction at the NASA 
Plum Brook Station. The test chamber, which is 100 feet in diameter and 122 feet high, 
will be surrounded by a 6-foot-thick concrete biological shield, parts of whose surface 
are to be borated. Of particular concern, from the materials activation standpoint, is 
the flux of thermal neutrons being returned by the concrete walls that resulted from 
incident fast neutrons from a test reactor in the chamber. 

Thermal neutron flux albedo factors resulting from fast neutrons are difficult to 
calculate, and there is a general lack of experimental data on the subject. Spinney 
(refs. 1 and 2) has used an age-diffusion theory approach to predict the albedo of a lightly 
capturing infinite-slab medium. This calculation tends to be unreliable since age- 
diffusion theory is questionable in boundary regions or media that have a noticeable 
hydrogen-atom density. Allen (ref. 3) and Leimdorfer (ref. 4) performed Monte Carlo 
calculations of neutron flux and current albedo factors for concrete slabs for source 
energies ranging over 0. 1 to 14 MeV. Because chemical binding effects and moderator- 
atom thermal motions are difficult to take into account and because machine time would 
become excessive, results of these calculations were cut off at 1 to 10 eV. Furthermore, 
only hydrogen, oxygen, silicon, and aluminum atoms are considered, and other elements 


in an assumed typical composition are replaced by silicon atoms. If the replaced atoms 
or others not known to be present cause appreciable absorption, particularly at low ener- 
gies, an accurate result cannot be predicted. 

Measurements of thermal -from -fast flux albedo factors resulting from fast neutrons 
were made for ordinary and borated concrete slabs by using a reactor beam at the Plum 
Brook Reactor Facility. Measurements were also made by using a collimated plutonium - 
beryllium source at the Lewis Research Center. A detector sweep method was used in 
both cases, foil detectors for the reactor source and a lithium 6 semiconductor sandwich 
detector for the plutonium -beryllium measurement. 


DEFINITION OF FLUX ALBEDO 


A distinction should be made as to the quantity of interest in consideration of albedo 
factors. In general the neutron albedo for a given material and geometry may be spec- 
ified in terms of current (or number), particle flux, energy flux, or dose. A neutron 
detector that is uniformly swept over a test slab surface has a response that is propor- 
tional to total emergent current, flux, or dose, depending upon detector characteristics. 
Thus foils that are thin (compared with the absorption mean free path) may be used to 
measure flux albedo and thick foils for current albedo. 

The quantity of interest is the thermal neutron flux albedo for a plane- collimated 
fast neutron beam incident upon an effectively infinite concrete slab. The thermal-from- 
fast neutron flux albedo |3rpy is defined as follows: 
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where an area, or surface, integration of emitted thermal (subcadmium) neutron flux 
cprjy over an area S is taken that corresponds to the surface integral of incident fast 
(epicadmium) neutron flux <p F> In practice, these surface integrals should be performed 
over a surface region large enough to account for most of the total emergent thermal 
neutrons. 
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Figure 1. - Experiment for reactor-beam measurement. 


EXPERIMENTAL PROCEDURE 
Reactor-Beam Source Experiment 

The first experiment on neutron flux albedo utilized the underwater beam room facil- 
ity at the NASA 60 megawatt Plum Brook reactor. The basic arrangement of equipment 
is shown in figure 1. A boral-filtered beam was incident upon the concrete slab surface 
over which gold and sulfur foils were swept. The concrete blocks that composed the test 
slab were poured according to the mix in table I. An approximate composition is given 
in table n, which is based upon the materials and proportions in table I. As is shown in 


3 



TABLE H. - NONBORATED A CONCRETE * 
COMPOSITION 

TABLE I. - MIX PROPERTIES FOR CONCRETE USED [Based ^ typical oxide analyses of 

IN REACTOR- BEAM EXPERIMENT cement and limestone aggregate. ] 


Constituent 

Composition, 
weight percent 

Calcium carbonate 

40. 6 

Silicon dioxide 

36.8 

Calcium oxide 

7. 3 

Water a 

6. 6 

Magnesium carbonate 

5.9 

Ferric oxide 

. 7 

Aluminum oxide 

1. 2 

Magnesium oxide 

.3 

Sulfur trioxide 

.3 

Sodium oxide 

. 1 

Potassium oxide 

. 1 

. Sulfur 

. 1 


a Water retention assumed to be 
100 percent. 


Material 

Nonborated A 

Borated A 


3 

Composition, lb/yd 

Cement 

(type II grey) 

470 

470 

Sand 

(Lake Erie) 

1357 

978 

Coarse aggregate 
(crushed limestone) 

2032 

2032 

Water 

272 

272 

Boron frit 

(Chicago Vitreous #S-85) 

... 

379 


figure 1, the gold and sulfur foils were mounted on a wand that was mechanically driven 
by two synchronous motors across an area approximately 18 inches square with an area 
resolution of about 0. 5 inch. Motor speed was such that a single sweep across the test 
area required a time of about 1 hour. Since the beam intensity was somewhat low, many 
sweeps (about 20) were required to obtain adequate sulfur pellet counting statistics. 

Also, since the beam intensity was subject to noticeable change over a reactor power 
cycle, all data were normalized to the activity of the sulfur pellet beam monitors placed 
in front of the boral shield (fig. 1). 

The necessary slab thickness required for saturation back diffusion of thermal neu- 
trons was determined by placing bare and cadmium- covered gold foils on the slab surface 
and by varying the number of concrete blocks in the slab. A typical plot is given in fig- 
ure 2 for the nonborated blocks, which indicates that a slab thickness of 12 inches or 
greater is effectively infinite for test simulation purposes. 

The thermal-from-fast flux albedo /3^,p is related to the foil activities by the 
following equation: 





( 2 ) 
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Figure 2. - Gold foil activity as function of slab thickness. 
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and 

A cadmium-difference gold foil specific activity 

5 

A„ sulfur pellet specific activity 

s 

e sulfur pellet counting efficiency 

s 

e gold foil counting efficiency 

t exposure time 

t^ decay time 

— 32 -4 2/ 

2 g S (n, p) macroscopic cross section averaged over E > 0. 5 eV, 4. 65x10 cm /g 

— 197 2 

S Au (n,y) macroscopic cross section averaged over 0 to 0. 5 eV, 0. 214 cm /g 

D 

X phosphorus 32 decay constant 

S 

\ gold 198 decay constant 

D 

The gold foil activity A^ was corrected for beam thermal neutron background by 

correcting for activity produced by sweeping foils with no concrete in the beam path. 

The mean sulfur cross section 2 0 corresponded to a microscopic value of 26 millibarns, 

s 32 

which was obtained by averaging S (n, p) cross-sectional data over a diffusion- code 
spectral- shape calculation for the beam (fig. 3) for neutron energies greater than 0. 5 eV. 
In order to obtain the mean gold activation cross section over 0 to 0. 5 eV, where 
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Figure 3. - Calculated beam spectrum. 

0. 5 eV is considered the cadmium cutoff, a Maxwellian density spectrum was assumed 
for 0 to 0. 1 eV, and a 1/E flux spectrum for 0. 1 to 0. 5 eV. The resulting mean cross 
section is 73 percent of the 2200-meter-per -second cross section value for a l/v de- 
tector such as gold. 

Measurements of for borated concrete were made by replacing the first layer 

(6 in. ) of nonborated blocks with borated blocks. These borated blocks contained approx- 
imately 0. 8 weight percent of elemental boron. 


Plutonium-Beryllium Source Experiment 

The experimental arrangement is indicated in figure 4, where a 5-curie plutonium- 
beryllium source was placed in a boral- shielded collimator to provide a fast neutron 
source. The concrete slab for this measurement was not composed of the same blocks 
used for the reactor source measurement, but consisted of limestone- aggregate struc- 
tural blocks. The same detector- sweep method of flux measurement was used as was 
used for the reactor-source measurement, except a lithium 6 semiconductor sandwich 
detector was used to count thermal neutrons. 

The total number of incident fast neutrons per unit time was taken as the product of 
the collimator exit area and the neutron flux and was determined by the source strength 
and an assumed l/r^ relation. Detector count rates were recorded with and without a 
20-mil cadmium cover. Cadmium-difference count rates were corrected for background 
by placing a large boral plate over the surface of the concrete slab. 

The flux albedo /3pp is given in the following equation: 
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where 

<p T surface- averaged thermal (E = 0 to 0. 5 eV) neutron flux 
surface area swept out by the detector 
cp-p fast (E > 0. 5 eV) neutron flux at collimator exit 
S c collimator exit area 

The relation of cadmium-difference detector count rate Crp to cp^ is given by 

Crp = KCTj^,j, (5) 

— fi 

where = Li (n, a) cross section averaged over 0 to 0. 5 eV and K is the detector 
constant. 

The flux albedo /3 TF may be written then in terms of measured parameters: 


/ 3 tf - 
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where Q is the source output and d is the distance from source to collimator exit. 
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RESULTS AND DISCUSSION 

The values obtained for albedo are summarized in table m. The value of /3 TF for 
the plutonium-beryllium measurement was 0. 13±0. 04. The indicated error includes 
counting statistics and uncertainties in both the lithium 6 cross section and mass. The 
concrete used in obtaining this value was a nonborated limestone aggregate, but was not 
the same concrete as that used for obtaining the reactor source value. The plutonium - 
beryllium albedo value includes an estimated correction factor to account for detector 
readings being recorded over a finite surface area about the beam center and thus not 
recording the total emergent thermal neutron flux. This correction is described in 
appendix A. 

Table HI also includes an analytical current albedo value of 0. 15±0. 10, which was 
calculated for ordinary concrete with an age-diffusion method (see appendix B). Since 
the experimental and analytical albedo values are on a flux or current basis, respectively, 
some assumption on emergent angular distribution of neutrons must be made for a valid 
comparison. The measured plutonium-beryllium flux albedo implies a current albedo 
of 0. 09±0. 03 if a cosine distribution for emergent thermal neutrons is assumed, since 
cos 6 is 2/3, where 6 is the angle from the outward normal to the surface. Thus the 
measured plutonium -beryllium albedo and the age-diffusion result are in agreement, 
within the ranges of error. 

The use of foil data for the reactor-beam source and equations (2) and (3) produce a 
value for /3 TF of 0. 020 for nonborated concrete. The cadmium-difference gold foil 


TABLE m. - NEUTRON ALBEDO VALUES 


Albedo 

Basis 

Source 

energy 

Concrete 

Neutron 

albedo 

value 

Flux 

Measurement 

Plutonium- 

Nonborated B 

0. 13±0. 04 



beryllium 



Current 

Measurement 

Plutonium- 

Nonborated B 

0. 09±0. 03 


estimate 3. 

beryllium 



Current 

Calculated* 5 

Fission 

Nonborated 

0. 15±0. 10 

Flux 

Measurement 

Reactor 

Nonborated A 

0.020 

Flux 

Measurement 

Reactor 

Borated A 

<0. 003 


Si 

Based upon measured flux value assuming a cosine emergent 
distribution. 
b See appendix B. 
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activity for the borated concrete, however, was essentially the same as that for no con- 
crete in the beam path. Since this beam thermal neutron background corresponded to a 
value of 0. 003 for 0^, the actual albedo for borated concrete would be smaller. No 
successful measurements of beam spectrum were possible since beam gamma intensity 
was too high to use lithium 6 semiconductor sandwich detectors, and beam neutron inten- 
sity was too low to use foils adequately. Consequently, a major source of uncertainty 
exists in the value for £ used in equation (3). 

The large difference in values obtained for the reactor and plutonium-beryllium 

sources is probably not accountable in terms of source spectrum effects; however, a 
large uncertainty in the reactor source albedo value arises from the lack of beam spec- 
trum data and the consequent error in mean sulfur cross section used for fast neutron 
dosimetry. Note that no flux-area correction factor (appendix A) was applied to the 
reactor source albedo value because of the lack of data. On the basis of the plutonium- 
beryllium correction, however, this effect would cause about a factor of 2 increase in 
the reactor source albedo value. The quoted reactor source albedo value could be in 
error by as much as a factor of 10 because of these combined uncertainties. 

The difference between albedo values for reactor and plutonium-beryllium sources 
might also be in part due to compositional variation in the concretes used for each meas- 
urement. The relative elemental thermal neutron absorption contributions, based upon 
the nonborated concrete A composition given in table II, are presented in table IV. Ap- 
proximately 93 percent of the total thermal neutron absorption is due to calcium, hydro- 
gen, and silicon. Calcium should also be a major absorber in the nonborated concrete 


TABLE IV. - RELATIVE THERMAL NEUTRON ABSORPTION FOR 
NONBORATED A CONCRETE 


Element 

Composition, 

weight 

percent 

Cross section, 
b 

Macroscopic 
cross section, 

S a> 

cm 2 /g 

Percent of 
total 2 

a 

Hydrogen 

0.7 

0. 332 

140 xlO 5 

37.4 - 

Oxygen 

51.0 

.0002 

00 

CO 

. 1 

Silicon 

16.2 

. 16 

56 

15.0 

Calcium 

23.0 

.44 

150 

40.2 

Aluminum 

.6 

.23 

3.1 

.8 

Iron 

. 5 

2. 53 

14 

3.7 

Magnesium 

1.9 

.063 

3.0 

.8 

Sodium 

. 1 

. 505 

1.3 

.3 

Potassium 

. 1 

2.07 

3.2 

.8 

Carbon 

! 5.7 

.0034 

1.4 

.4 

Sulfur 

.2 

. 52 

2.0 

. 5 
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B, used for plutonium-beryllium measurements, since this was also a limestone- 
aggregate concrete. 


« 


SUMMARY OF RESULTS 

Measurements of thermal-from-fast neutron flux albedo were made for non- 

borated and borated concrete by using a reactor-beam source. Flux albedo measurements 
were also made on a nonborated concrete with a collimated plutonium-beryllium source. 
The following results were obtained: 

1. On the basis of reactor-beam measurements, a thickness of 12 inches of nonbo- 
rated limestone concrete was found to be sufficient for essential saturation back diffu- 
sion of thermal neutrons. 

2. The best value obtained for /3^,p is 0. 13±0. 04 for plutonium-beryllium neutrons 
incident on a nonborated concrete. This value is in agreement with an age-diffusion 
theory calculation for ordinary concrete. 

3. The reactor-beam experiment yielded a value of 0. 020 for from nonborated 

concrete. A value of /3rpp for borated concrete was below sensitivity limits, which 
corresponded to an albedo value of 0. 003. The large difference between quoted plutonium- 
beryllium and reactor source albedo values might result from a large uncertainty in 
reactor-beam source spectrum, incomplete flux-area integration, and composition var- 
iation that may exist for the two different nonborated concretes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1966, 

120-27-01-03-22. 
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APPENDIX A 


FLUX-AREA CORRECTION FACTOR FOR PLUTONIUM-BERYLLIUM 

SOURCE MEASUREMENT 


The lithium 6 detector thermal neutron readings were recorded over the dashed-line 
square areas indicated in figure 5, where d is 18 inches. To account for thermal neutron 
flux beyond the areas of actual measurement, a simple analytical model was assumed 


(fig. 6), where the plutonium-beryllium neutrons 




Figure 6. - Point-source model. 


incident upon the concrete were assumed 
to generate a unit point source of ther- 
mal neutrons at a depth equal to an 
average slowing-down length L g . The 
thermal neutron flux at the surface is 

-Kp 

<p( r ) = — — — (Al) 

4ttp 2 

where the coefficient K is evaluated 
from data recorded over finite surface 
areas. 

The relative surface flux-area 
integrals enclosed by radii r ^ and rg 
are first determined by a surface geom- 
etry transformation between square 
areas used for measurements and cir- 
cular areas defined by the model. The 
values of r^ and are determined 
by equating the area d 2 to the central 
circular area and also to that of the 
segment defined by 6 (fig. 5). If a 
is defined as the ratio of flux-area 
integrals as measured for the area of 
squares B and A (fig. 5), then a * , 
which represents the ratio of flux-area 
integrals for the annular region of 
dimension r 2 - r i to that for the 
region inside r^, is 
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a' = *la 


(A2) 


where 6 equals 2 tan"* [d/( r i + r 2 )] an< * * s determined from the geometry transformation. 
A value for K is determined by solving the following equations: 


f 


<p( r)r dr 


f r \ 

I <p(r)rdr 
* / 0 


_ E l( K V r l + L s) - E l( K V 

Ei(Kt s ) - Ej( K ^f 


2 t2 
r 2 + L s 




(A3) 


A flux-area correction factor y is defined as the relative surface integral of flux 
outside dimension r 0 : 


E l( K V 

y = 1 L_ 


2 t-2 

r 2 + L s 


E 1 (KL g ) 


(A4) 


The corrected value of flux albedo 0, which includes flux-area contributions from 
all surface area regions, may be expressed as 


0 = 0Q + + VP 

= !+“! a. 


(A5) 


where 0 q is the value of flux albedo based upon measurements over the area of the center 
square A (fig. 5). The quantity a' is experimentally evaluated and fitted to the model 
while y is evaluated with the model. 

The following is a summary of parameters including the resultant value of 0 (the 
value of L g was assumed on the basis of representative values of age r from references 
2 and 5: 
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r i = 10. 2 in. 


rg = 29. 4 in. 


0 = 0. 855 rad 


L g = 5±2 in. 


a = 0. 082± 0.027 


a' = 0. 60±0. 20 


A A7 , 0. 13 

y = °- 07± o. 02 


| 3 0 = 0. 0764±0. 0184 
/3 = 0. 13±0. 04 

A disk source would have been a more reasonable source geometry to select for this 
model, but flux-area integral calculations for such a source cannot be evaluated simply. 
The validity of the point source model, however, can be investigated by displacing the 
point source a distance m (fig. 6) from the center of symmetry and by calculating the 
relative change in flux-area product over a surface area with a radius of 10 inches. This 
area is that included by dimension r^ (see fig. 5) and over which the value of /3 q is 
determined. The disk source was estimated from the experiment geometry to be about 
9 inches in diameter. The appropriate value of m is not the source radius, but corre- 
sponds to a radial distance averaged over the source area. Thus m was taken as two- 
thirds of the source radius, or 3 inches. It was estimated from the results of Rockwell 
(ref. 6) that less than a 5-percent change in flux-area product over the 20-inch-diameter 
surface area is effected by such a point-source displacement. Hence, the original point- 
source assumption is a reasonably valid source geometry for the surface area sizes upon 
which measurements were based. 
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APPENDIX B 




f 


ANALYTICAL ALBEDO CALCULATION 


An estimate of thermal-from-fast neutron current albedo (3j was made by using 
the diffusion theory result (ref. 2): 



SL 


(L g + L)^l + 


2D 

L 


(Bl) 


where 

S total source strength 

L thermal diffusion length in concrete 

L slowing- down length in medium appropriate to thermal neutrons 

D thermal diffusion coefficient 

The total source strength S is taken as 1 - F(t t1i ), where F(t t1i ) represents the 
fraction of neutrons emerging form the concrete with age less than thermal. 

An expression for F(r,j, h ) has been obtained (ref. 1) by using age theory with a 
distributed source of fast neutrons given by first collisions of incident beam neutrons. 

An estimate of F(t t1i ) was made by utilizing a set of curves (ref. 2) plotted as functions 
of the parameters r Th , and £ Tr , where is the thermal neutron age in con- 

crete, S R is the macroscopic removal cross section for fast neutrons incident upon 
concrete, and £ Tr is the transport mean free path for neutrons slowing down to thermal 
energies in concrete. 

The following is a summary of mean parameter values that corresponded to several 
ordinary concretes (ref. 5) (the parameter values assume a water retention of 100 per- 
cent, and parameter variation about the mean value was assigned on a confidence level 
of 90 percent): 


L = 8. 3±4. 2 cm 
D = 0. 66±0. 30 cm 
T Th = 187±98 cm^ 
L = 13. 7±3. 6 cm 

o 
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Jfrj, r = 1. 98±0. 60 cm 
S R = 0. 081±0. 005 cm -1 

F ( T Th> = °* 55±0 - 10 

When these values are used in equation (Bl), a value of 0. 15±0. 10 is obtained for /3j. 
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